Entanglement in Quantum Spin Chains, Symmetry Classes of Random Matrices, and 

Conformal Field Theory 
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We compute the entropy of entanglement between the first N spins and the rest of the system in 
the ground states of a general class of quantum spin-chains. We show that under certain conditions 
the entropy can be expressed in terms of averages over ensembles of random matrices. These averages 
can be evaluated, allowing us to prove that at critical points the entropy grows like « log 2 N + k 
as TV — > oo, where k and k are determined explicitly. In an important class of systems, k is equal 
to one-third of the central charge of an associated Virasoro algebra. Our expression for k therefore 
provides an explicit formula for the central charge. 
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Entanglement has recently come to be viewed as an im- 
portant physical resource for manipulating quantum in- 
formation. The problem of quantifying it is, however, still 
poorly understood, especially when the entanglement is 
shared between more than two systems. When the en- 
tanglement of a pure state is shared between two parties, 
i.e. in a bipartite system, Bennett et al [l| have shown that 
it is consistent to define it as the von Neumann entropy 
of either of the two parts. We consider here the gen- 
eral class of quantum spin chains arising from quadratic 
chains of fcrmionic operators in their ground state. These 
systems are partitioned into two contiguous subchains. If 
the ground state is non-degenerate, this subdivision cre- 
ates a pure bipartite system; our main result is to calcu- 
late its entanglement entropy by relating the problem to 
one in random matrix theory. 

As is well known, the systems we are studying exhibit 
quantum phase transitions. These manifest themselves 
as qualitative changes in the decay of correlations: alge- 
braic at a critical point and exponential decay away from 
it. Entanglement plays a fundamental role in the quan- 
tum phase transitions that occur in interacting lattice 
systems at zero temperature [E IE El IE IE • Under 
these conditions the system is in the ground state, which 
is also a pure state, and any correlations must be a con- 
sequence of the fact the ground state is entangled. It 
follows immediately that the entanglement changes qual- 
itatively at critical points. 

In this context, Vidal et al. |j| studied the ground 
states of a wide range of one-dimensional spin models 
partitioned into two consecutive subchains. They ob- 
served numerically that, when the Hamiltonian under- 
goes a phase transition, the entanglement of formation 
of these bipartite systems grows logarithmically with the 
size N of one of the two parts. Jin and Korepin then 
proved that the entropy grows like | log 2 N in the XX 
model, for which the Hamiltonian is 
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where <j a denotes the Pauli matrices and a = x, y, z. Re- 
cently, Korepin [5j and Calabrese and Cardy [(j showed, 
using conformal field-theoretic arguments developed by 
Holzhcy et al. |9j, that the logarithmic divergence of 
the entanglement in one dimensional systems is a gen- 
eral consequence of the logarithmic growth of the entropy 
with the size of the system at phase transitions. These ar- 
guments determine the constant multiplying the leading 
order log 2 -/V term in the asymptotics to be one-third of 
the central charge of the associated Virasoro algebra • 

We show here that if a quantum spin-chain Hamil- 
tonian posses certain symmetries, the entanglement can 
be expressed as an average over an ensemble of random 
matrices corresponding to one of the classical compact 
groups equipped with Haar measure, i.e. one of the fol- 
lowing groups: U(7V), Sp(2/V) and ± (N), where the 
superscript ± indicates the connected component of the 
orthogonal group with determinant ±1. From the point 
of view of the entanglement entropy (and of spin-spin 
correlations), quantum spin chains therefore divide into 
symmetry classes related to the classical compact groups. 
The XX model turns out to be an example of a system 
with U(iV) symmetry. The averages that occur can be ex- 
pressed either as Toeplitz determinants (i.e. determinants 
of matrices in which the elements are functions of the dif- 
ference between the row and column indices) , in the case 
of U(JV), or as determinants of specific combinations of 
Toeplitz and Hankel matrices (i.e. matrices in which the 
elements are functions of the sum of the row and col- 
umn indices) for the other compact groups. Asymptotic 
formulae for these determinants then lead to general ex- 
pressions for the leading-order and next-to-leading-order 
terms in the asymptotics of the entanglement in the limit 
as the total number of spins tends to infinity and then as 
N — > oo. 

We find that at a critical point the entanglement grows 
logarithmically with N, in agreement with the conformal- 
field-theoretic calculations. We derive a general formula 
for the associated constant of proportionality. This is 
a rational number, the numerator of which is shown to 
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factorize into a universal part, related to symmetries of 
the quantum Hamiltonian and which can be calculated 
from the random-matrix averages, and a non-universal 
(i.e. Hamiltonian-specific) part, which we also evaluate. 
In the unitary case, comparing with the results in 
leads to an explicit formula for the central charge. How- 
ever, our approach also extends to systems where the 
conformal-field-theoretic results cannot be applied di- 
rectly. These are the systems related to the other com- 
pact groups. Further details of our calculations and re- 
sults may be found in [ll| . 

The most general form of Hamiltonian related to quan- 
tum spin chains is 



block structure 



H a = a 
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where a and 7 are real parameters, < 7 < 1, the bjS are 
Fermi oscillators, A is an Hermitian matrix, and B is an 
antisymmetric matrix. We take periodic boundary con- 
ditions, i.e. bu = bo- Without loss of generality, we will 
consider only matrices A and B with real elements. The 
Hamiltonian @ can always be re-expressed in terms of 
the Pauli spin matrices using the Jordan- Wigner trans- 
formation [U03- 

We will here be concerned with the entanglement be- 
tween the first N oscillators and the rest of the chain 
when the system is in the ground state |* g ) and as the 
length of chain tends to infinity. We decompose the 
Hilbert space into the direct product H = Hp <£> Hq, 
where Hp is generated by the first N sequential oscilla- 
tors and Hq by the remaining M — N. Our goal is to 
determine the asymptotic behaviour for large N << M 
of the von Neumann entropy Ep = — Tr pp log 2 pp , where 
Pp = Tr Q p PQ and p PQ = |* g )(* g |. 

The first step involves determining the expectation val- 
ues with respect to |^ g ) of products of arbitrary numbers 
of the operators 
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From the invariance of the Hamiltonian J3J) under the 

:* g |™;|* g ) = 



transformation 



—bj, it follows that 



0; for the same reason, the expectation value of the 
product of an odd number of m^s must be zero. The 
expectation values (*& g \mjmk can be deduced us- 
ing the approach of Lieb et al [13]: (^ g \m,jmk |* g ) = 
8jk + i(^Af)jfc) where the correlation matrix Cm has the 
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with 
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the matrix Tm is defined by 

M-l 

( T M) jk = ^ik, j,k = 0,...,M-l, (6) 
1=0 

and the vectors <f> k and xp k are real and orthogonal and 
obey the eigenvalue equations 

a 2 U-h- 7 b\ (a - h + 7 b\ 4> k = |A fe | 2 fe , 



a 2 (A-|/4 -/;)(. I --/ -:B)th, =|A,j- ( /',. 
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These vectors are related by 



a[A--I + jB)ct> k = \A k \il> k , 



a{A--I- 1 B)xP k = \A k \cf> k . 



(8a) 
(8b) 



The expectation values of the product of an even number 
of rrijS can then be computed using Wick's theorem (see, 
for example, 0). 

Following the calculation in 0, the formula for the 
entropy of the subchain P that one obtains using these 
expressions for the expectation values is then 



Ep = lim lim 

e^o+ <5-»o+ 2-7ri 



where 



c(e,S) 



e(l + e, A) ; " v ' dA, 



, s x + v ( x + v 
e{x, v) = — log 2 



dA 



x — V I x — V 



(9) 



(10) 

-Ojv(A) = det (IX — S), S is the real symmetric matrix 

1/2 _ 
(TprJ%) and Tjv is obtained from the matrix JBJ by 

removing the last M — N rows and columns. The contour 

of integration c(e, 5) depends on the parameters e and 

8 and includes the interval [—1,1]; as e and 8 tend to 

zero the contour approaches the interval [—1,1]. This 

guarantees that the branch points of e(l+e, A) lie outside 

the contour of integration and thus that e(l + e, A) is 
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analytic inside c(e, S). The eigenvalues of S must all lie 
in the interval [—1,1], and this is the case for the various 
Hamiltonians we consider [ll|. 

We first specialize to cases where the Hamiltonian 
is invariant under translations. For example, the XX 
model has this symmetry. 

We denote A = a A — 21 and B = a^B. If H a is in- 
variant under translations of the lattice {0, 1, ... , M— 1}, 
then the elements of the matrices A and B must depend 
only on the difference between the row and column in- 
dices, i.e. A and B must be Toeplitz matrices. In addi- 
tion, because of the periodic boundary conditions, A and 
B must be cyclic. 

Now, let a and b be two real functions on Z/MZ, even 
and odd respectively. The matrix elements of A and B 
can be written as 

Ajk = a(j - k) and B jk = b(j - k). (11) 

The complex exponentials 



M 



j,k = Q,. 



,M-1, 



(12) 



form a complete orthonormal set of eigenvectors of cyclic 
matrices, as can be easily verified by direct substitution. 
The matrices A and B defined in (|11J) commute. As a 
consequence, the complex exponentials <|12|) are a com- 
plete set of eigenvectors of A + B too. 

The eigenvalues of A + B can be determined by insert- 
ing the eigenvectors <(T2)l into the eigenvalue equations 
and using the parities of the functions a(j) and b(j). We 
have that when M is odd 

(M-l)/2 

A fc =a(0) + 2 Hj) cos kj + ib(j) sin kj] (13) 

and when M is even 



A fc = o(0) + (-l)'a(M/2) 

M/2-1 

+ 2 51 [a(j) cos fcj + sin fcj] , 

3=1 



(14) 



where fc = 2wl/M. Then 
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where A(0) is the periodic function 

CO 

a(0)= y, 



-Ki-W<i9, (15) 



(16) 



with Aj = a(j)—b(j) if j > and Aj = a(j)+b(j) if j < 0. 
Note that T N [g] is a Toeplitz matrix. g{9) = A(O)/\A{0)\ 
is called the symbol of T. It is worth emphasizing 



that (|15|) has been obtained by assuming only the trans- 
lation invariance of the Hamiltonian (|5J and periodic 
boundary conditions. 

We now make the key observation that when 7jv[g] 
is symmetric the matrix S, which appears in the defi- 
nition of Dn(X), is equal to Tm- We can then apply a 
famous identity of Heine [T^J and Szego 0], which as- 
serts that if G(U) is a function on V(N) that depends 
only on the eigenvalues exp(i8j) of U and is such that 
G{U) — YljLi 9$ j)i where g(8) is 27r-periodic, then 



G{U) 



U(JV) 
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j- k >j,k=0,...,N-l ' 
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where gi is the ith Fourier coefficient of g. In our context, 
this implies that © can be expressed as an average with 
respect to Haar measure over the unitary group U(A r ), 
i.e. over the Circular Unitary Ensemble of N x N random 
matrices. A necessary and sufficient condition for T^[g] 
to be symmetric is that A(6) should be real and even, or 
equivalently 7 should be zero; in other words, the inter- 
action in the Hamiltonian (J2J must be isotropic. When 
7 = the symbol g(9) is a piece- wise continuous function 
which takes the values 1 and —1 and has discontinuities 
at all points 8 r where the equation 



A(9 r ) = 



(18) 



is satisfied, with the additional condition that the first 
non-zero derivative of A(9) at 9 r is odd. 

Given that under the general conditions specified 
above the entropy of entanglement can be expressed as 
an average over U(AT), it is natural to ask whether under 
different conditions it can be expressed as an average over 
random matrices drawn from the other classical groups. 
The question is: how are the symmetries of the Hamilto- 
nian reflected in the group which determines the entropy 
of entanglement? 

We begin with the orthogonal group + (2N). The 
analogue of the Heine-Szego identity in this case relates 
group averages to the determinant of a sum of Toeplitz 
and Hankel matrices. A straightforward calculation gen- 
eralizing that given above shows that this can be ar- 
ranged for the spectral determinant Djy if (and only if) 
7 = and Aji — a(j — I) + a(j + I) where, because of 
the periodic boundary conditions, a must be a function 
on Z/MZ and must also be even in order for A to be 
symmetric. Note that the Hamiltonians in this class are 
not translation invariant. The properties of Dn are the 
same as those in the unitary case, except that T?{[g] is 
the sum of a Toeplitz and a Hankel matrix. The symbol 
has the same general form as in the unitary case. 

The calculations for the other compact groups follow 
exactly the same pattern except that for Sp(2A) and 
0-(2N + 2),_A jk = a(j - k) - a(j + k + 2), and for 
± (2iV + l), A jk = a(j-k)Ta(j + k + l). A gain, in these 
cases the Hamiltonians are not translation invariant. 



4 



The asymptotics of the entropy of entanglement can 
now be calculated using the Fisher-Hartwig conjecture 
for the determinant Djy in the unitary case UM and 
recent generalizations of this conjecture in the other 
cases |l7l [lg , and then by computing the integral in J3J . 
The result is that as N — > oo 



2 W °R 
6 



log 2 N, 



(19) 



where R is the number of solutions of (|18fl in the interval 
[0, 7r) and 

{1 if the average is over U(iV) 
otherwise. 

The asymptotic relation (|19|l represents our main re- 
sult. In the unitary case, comparing with the results 
of nsn, it provides an explicit formula for the central 
charge, which may be seen to depend in a non-trivial way 
on the geometry of the Hamiltonian. In the case of the 
other classical compact groups, when the Hamiltonian 
is not translation invariant, the conformal- field-theoretic 
results do not apply directly. The factor 2 WG is univer- 
sal, depending only on the symmetries determining the 
classical compact group to be averaged over. The factor 
R is Hamiltonian-dependent. For the XX model, which 
is an example with unitary symmetry, R = 1 and l|19|) 
coincides with the formula derived in |4J. 

Lower order terms in the Fisher-Hartwig conjecture 
and its generalizations lead directly to general formulae 
for the next-to-leading-order (constant) term k in the 
asymptotics of the entropy of entanglement when N — > 
oo. For the unitary group we find 



KU(iV) 



R 



3 In 2 



(K - 6/ 3 In 2) 



(21) 



where -Z3 is a constant evaluated in [4| to be 0.0221603 . 
and 



X = l + 7E + i^ln|l-e i2 ^ 



I E (-D (r+S) ln 



Kr<s<R 



(22) 



Here 7e is Euler's constant. When R = 1 this equa- 
tion reduces to the result of Jin and Korcpin for the XX 
model. For the other compact groups we hnd, similarly, 
that R — ku(tv)/2 + R/Q. 

We end with some general remarks. First, the circu- 
lar ensembles of random matrices may be seen to play a 
special role in the context of the spin chains and bound- 
ary conditions we have considered here. It would be 
interesting to know whether the other random matrix 
ensembles may be used to describe systems with differ- 
ent interactions and boundary conditions. Second, we 



note that away from critical points A(8) is continuous, 
and so R = 0, which is consistent with previous observa- 
tions that the logarithmic growth of Ep is a critical phe- 
nomenon. Our approach determines the limiting value 
of Ep away from critical points too. Finally, the calcu- 
lations and results described above extend straightfor- 
wardly to spin-spin correlations in the families of quan- 
tum spin chains we have considered. 

We gratefully acknowledge discussions with Estelle Ba- 
sor, Peter Forrester and Noah Linden, and the kind hos- 
pitality of the Isaac Newton Institute, Cambridge. 
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